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^! We study the relevance of strange degrees of freedom for nucleon structure functions. 

Oh: 



For this purpose we employ the three flavor generalization of the collective quanti- 
zation approach to the chiral soliton of the bosonized Nambu-Jona-Lasinio model. 
Contrary to many other soliton models the hadronic tensor is tractable in this model. 
By applying the Bjorken limit to the hadronic tensor we extract the leading twist 
contributions to the nucleon structure functions at the low energy scale at which the 
model is assumed to approximate QCD. After transforming to the inflnite momentum 
frame and performing the DGLAP evolution program to these structure functions we 
compare with available data for deep inelastic electron-nucleon scattering. 
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1. Introduction 

In this paper we present the computation of the leading twist pieces of both the unpolarized 
and polarized nucleon structure functions in the three flavor Nambu-Jona-Lasino (NJL) [Q] 
chiral soliton model 0, 0, Q. Applying functional bosonization techniques this model for 
the quark flavor dynamics can be straightforwardly expressed in terms of meson degrees of 
freedom. The identification of baryons as solitons of effective meson theories is motivated 
by generalizing quantum chromo dynamics (QCD) to a large number {Nc) of color degrees 
of freedom |^. Essentially this generalization shows that the baryon properties exhibit the 
same dependence on A''^ as expected for soliton solutions of the effective meson theory. In 
addition to this analogy there are various reasons based on phenomenological observations to 
consider soliton models suitable for studying baryon properties. For example, such models quite 
successfully describe many aspects of static baryon properties [^ P, ^ [l^ . This is in particular 
the case for observables related to the proton spin puzzle |jll| which is directly related to the 
polarized structure functions of the nucleon. Working in the three flavor extension of the model 
provides us with the advantageous feature of direct access to the strange quark contributions 
to structure functions even though the net strangeness of the nucleon vanishes. Hence we can 
gauge the importance of these degrees of freedom for dynamical properties of the nucleon. 
This is particularly interesting because the analyzes of experiments performed at CERN and 
SLAC have indicated that the spin of the strange quarks inside the nucleon yields a sizable 
contribution to the total nucleon spin. 

Unfortunately there are a number of soliton models available. The reason for adopting the 
NJL chiral soliton model lies in the fact that the main ingredient in any structure function 
calculation is a bilocal commutator of hadronic currents as can be observed from the defining 
equation ( p4|) of the hadronic tensor. It appears that its calculation in models with elementary 
meson fields is infeasible. However, it is much less involved in a model defined in terms of 
quark fields because the Compton amplitude, whose absorptive part equals this commutator, 
can be traced in the process of functional bosonization. Formally this process provides a unique 
regularization of the quark loop. Leaving aside the cumbersome issue of regularization in the 
NJL model, the leading twist contribution to the absorptive part of the Compton amplitude 
becomes as simple as a bilocal and bilinear operator of the quark fields. 

In the present approach to the three flavor NJL-model we do not include the 't Hooft de- 
terminant [0 which, when included [0, provides a mass for the singlet component of the 



ry-meson and parameterizes the mixing of the uncharged pseudoscalar mesons. The reason for 
ignoring this piece of the action is that pseudoscalar flavor singlets couple only weakly to the 



chiral soliton, if at all |Tl]]. In the soliton approach the influence of the 't Hooft determinant on 
nucleon properties is hence negligible. It should however be remarked that at a sub-leading or- 
der in 1/Nc this term might give some contributions of the strange quarks to nucleon properties 
due to a mixing of non-strange and strange constituent quarks [Q. 

We should mention that there is an exhaustive number of model calculations for nucleon 



structure functions. The models applied range from various kinds of bag models |TB|, [T^ over 
cloudy bag models [0 and meson cloud models ITsf to (diquark) spectator models [|l^, just 



to name a few. In the context of structure functions those model calculations may eventually 
be somewhat more reliable than soliton model results. After all, many of those models have 
particularly been set up to study structure functions. On the other hand soliton models are very 
attractive not only because of the above mentioned connection to QCD in the limit of large Nc 
and the possibility to generalize them to three flavors but over and above they provide a very 
comprehensive picture of baryons ranging way above static properties of nucleons from meson- 
nucleon scattering |2^, ^ over systems with a heavy quark ||2l| to applications in the context 



of nuclear matter |2^. Thereby these models use only very a limited number of parameters. 



The present study apparently serves to further complete this comprehensive picture. 

This paper is organized as follows: In section 2 we briefly review the NJL chiral soliton model 
and discuss the inclusion of strange degrees of freedom within the collective approach. Effects 
due to flavor symmetry breaking are incorporated into the collective nucleon wave-function by 
a generalization of the Yabu-Ando method [^, [10|. We then explain the calculation of the 



hadronic tensor for localized field configurations [15]. Employing pertinent projection tech- 



niques the structure functions are extracted from this tensor. In this section we furthermore 



include a brief discussion on the relevance of boosting to the infinite momentum frame p4| , p5[| . 
In the 1/Nc expansion scheme the dynamical response of the soliton to the momentum trans- 
fer is omitted which results in a violation of Lorentz covariance. Fortunately these omissions 
become irrelevant in the infinite momentum frame. We also discuss the valence quark approxi- 
mation ^^ to the structure function calculation which is to be considered as an approximation 



to avoid ambiguities inevitably arising when regularizing the vacuum part to bilocal-bilinear 
quark operators. Some of the previous approaches |2^, ^ to the two flavor structure functions 
adopted a once-subtracted Pauli-Villars type scheme. This, however, is not sufficient to render 
finite the quadratic divergences contained in the model. As a consistent incorporation of fiavor 
symmetry breaking in the three fiavor model demands the solution of the quadratically diver- 
gent gap-equation a single subtraction will not be suited for present studyS. In section 3 we 
present our results for the unpolarized structure functions for both the electric and the strange 
current. Subsequently we discuss the polarized structure function gi in section 4. In particular 
we refiect on the electric and strange pieces in the light of the proton spin puzzle. In section 5 
we will perform the DGLAP program to evolve the predicted structure functions from the low 
energy scale at which the model supposedly approximates QCD to the energy scale encoun- 
tered in the respective experiments. In section 6 the numerical results will be presented and 
compared to data from CERN and SLAG experiments. Finally section 7 serves to summarize 
and conclude our study. Technicalities of the calculations are relegated to appendices. 

Some of the results presented in this paper on the strangeness contribution to the polarized 



structure function g\{x) have already been published earlier ||29 



^In the case of the structure functions the absorptive part of the Compton amplitude - a matrix element 
of an operator which is quartic in the quark fields - appears to be the suitable starting point for a consistent 
regularization with two subtractions. 



2. The Hadronic Tensor from The Chiral Soliton in the 
Three Flavor NJL Model 

We split up this section into two pieces. In the first one we will review the treatment of the 
chiral soliton in the three flavor NJL model while the second one contains the discussion of the 
hadronic tensor Wfj_i, in this model. 

2a. The Chiral Soliton in the Three Flavor NJL Model 

The basis of our considerations is the NJL model Lagrangian [|l| 

C = q{z<^ - m')q + 2Gnjl E ( (^y?)' + (^^^75?)'! (1) 

which is chirally symmetric for a vanishing current quark mass matrix m°. In eq (^ q denotes 
the quark field. The matrices X'^/2 are the generators of the flavor group U{Nf). We will 
consider the case of three flavors {Nf = 3) and neglect isospin breaking, i.e. the current quark 
mass matrix acquires the form m° = diag(m°, m^, m^) as r7i° = m° = m^ is assumed. The NJL 
couphng constant Gnjl will later be determined from meson properties. 

Applying functional integral bosonization techniques the model (P can be rewritten in 
terms of composite meson fields [|]. The corresponding effective action is given by the sum 
•^NJL = -^F + -^m of a fermion determinant 

Af = Tr \og{ip) = Tr log (i^ - {PrM + PiM^f)] (2) 

and a purely mesonic part 

A^= /c/^xf-— ^tr(MtM-m°(M + Mt) + (m°)2)y (3) 

J \ 4Gnjl / 

Here Pr^l = (1 ±75)/2 denote the projectors onto right- and left-handed quark fields, respec- 
tively. The complex matrix M = S + iP parameterizes the scalar and pseudoscalar meson 
fields. 

Obviously, the action Af is quadratically divergent and needs to be regularized. We will 
apply the 0(4) invariant proper time regularization to Af in Euclidean space which then 
becomes a complex quantity Af = Ar + Aj. Proper time regularization amounts to replacing 
the real part of the fermion determinant Ar by a parameter integral 

Ar = ^Tr log (4lfe) ^-l Q^ y Tr exp (-s^^i.) (4) 

which becomes exact in the limit A — ;> oo up to an irrelevant additive constant. The imaginary 
part 

-A, = ^Trlog((4)^'^i.) (5) 



is finite and does not undergo regularization in order not to spoil the anomaly. 

To study the three flavor extension of the model it is suitable to introduce the following 
parameterization for the meson fields 

M = ^S^ (6) 

The matrix S is Hermitian whereas ^ is unitary. In this parameterization the chiral 811^(3) x 
SU/j(3) transformation M -^ LMK^ with constant matrices L and R is realized via S -^ KUK'^ 
and ^ -^ L^K^ = KC,R^ . The latter equation determines the unitary matrix K ^0 . 



Varying the action with respect to S yields the Schwinger-Dyson or gap equation which 
determines the vacuum expectation value (VEV), (S) = diag(m, m,ms) and relates it to the 
quark condensate {qq)i 



rrii 



,0 , ™3^c;Gnjl / ,rr^-^2\ _ _o nr^ It,. 
27r2 I ' ^ A 



< + <^^^^r?^r -1, C-^r =m^- 2Gnjl gg , (7) 



Here rrii = m, rUs denote the constituent masses of non-strange and strange quarks, respectively. 
The appearance of an incomplete Gamma-function of the order —1 reflects the fact that the 
model is quadratically divergent in the ultraviolet. 

Space-time dependent fluctuating pseudoscalar meson fields rjai^x) are introduced via 



^/(x) = exp ( z ^ ria{x)\a/2 j . 



Expanding the effective action up to second order in the fluctuations r]a{x) allows us (after 
rotating back to Minkowski space) to extract the inverse propagator for the pseudoscalar mesons 

Dij^kM ) = ( 76^ ^^M )%ihr (9) 

The polarization operator njj(g^) is given by 

n,(g^) = -2g^/g.(g^) + 2K - ^,ffU<i') - ^K - -?)(^ - ^) (lo) 

wherein 

fW) = \im^ + ^3?^ [ dx r(0, [(1 - x)ml + xm] - x(l - x)q^]IX^) . (11) 

The Bethe-Salpeter equation D^^{q^)P = which determines the physical meson masses mphys 
is equivalent to the condition that the meson propagator acquires a pole: 

det(Dr.\,(g2 = mJ,yj)=0. (12) 

Note that ffjiq"^ = ^phys) is the corresponding on-shell meson decay constant. From eq ([Tl| ) 
we can read off the pion decay constant 

f, = m^^ /' dxT{0, [m' - x(l - x)ml]/A') (13) 



47r2 
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Table 1: The up and strange quark constituent and current masses, the cutoff and the pion and 
the kaon decay constants for the parameters used later. Results are taken from ref pT| . 



m (MeV) 


m, (MeV) 


mVml 


A (MeV) 


U (MeV) 


Ik (MeV) 


350 


577 


23.5 


641 


93.0 


104.4 


400 


613 


22.8 


631 


93.0 


100.3 


450 


650 


22.4 


633 


93.0 


97.4 


350 


575 


24.3 


698 


99.3 


113.0 


400 


610 


23.9 


707 


103.0 


113.0 


450 


647 


23.6 


719 


105.7 


113.0 



as well as the kaon decay constant 



fl 



-{m + m. 



' 47r2 



dxTyO, [xni^ + (1 — x)rn?g — x{l — x)m|^]/A^j. 



(14) 



Using eqs (|12|,|13|) the four parameters of the model (|1|), the coupling constant Gnjl, the cutoff 
A and the two current masses mP and m° may be determined. The current quark masses are 
obtained from the pion and kaon masses, m^ = 135 MeV and itlk = 495 MeV. Fixing now the 
pion decay constant /^r = 93 MeV yields too small a value for the kaon decay constant, see 
table |I|. On the other hand, requiring fK = 113 MeV leaves us with too large a value for /^.H 
In the following we will exclusively employ parameter sets which reproduce the physical value 

of /,r. 

The expressions for the decay constants are only logarithmically divergent, which in a Pauli- 
Villars scheme could be cured by a single subtraction. For the full three flavor model this, 
however, would not be sufficient because the pertinent treatment of flavor symmetry breaking 
also requires the solution to the quadratically divergent gap-equation (|^. In ref ||32| the sizable 
effects of the second subtraction are discussed for meson observables which can be made finite 
by a single subtraction. 

For the investigation of the baryon sector we constrain the meson fields to the chiral circle, 
i.e. we replace S by its VEV (S) in eq (^). Hence our NJL model soliton has the more 
complicated structure 

M{x)=US){t)Ux) (15) 



with 



U^) 



exp [iQ{r)f ■ r/2] 




(16) 



containing the hedgehog ansatz. The computation of the classical energy associated with the 



^Determining /^ or Jk fixes the ratio K/m. This leaves one adjustable parameter, e.g. the coupling constant 
Gnjl- However, as Gnjl is not very transparent we will use the gap equations (|7|) to re-express it in terms of 
the up constituent mass m. 



unit baryon number sector of this field configuration is standard |3^, ^ 

^[9] - E[e = 0] = ^6v (1 + sign(6v)) + ^ 1^ ^J^ E e^P h^l 

+m2/2|rfV[l-cose(r)], (17) 

where e^ are the eigenvalues of the single particle Dirac Hamiltonian 

ho = a ■ p + m l3 expU'y^T ■ f<d{r))T + m^ (3S (18) 



In particular ey refers to the valence quark level which is strongly bound in background of the 
hedgehog. From the appearance of the strange and non-strange projectors S = diagf)(0,0, 1) 
and T = diagfl(l, 1,0) we observe that the strange quarks are not effected by the hedgehog field 
and hence obey a free Dirac equation. The profile of the chiral angle 0(r) is determined by 
self-consistently minimizing the energy functional ( p!7D 0, 3i . 



In order to generate states with baryon quantum numbers we apply the collective approach. 
In this approach time dependent coordinates A{t) G SU{3) which parameterize the orientation 
of the soliton in flavor space are introduced to approximate the unknown time dependent 
solution to the full equations of motion, i.e. 

ax,t) = A{t)Ux)AHt) ■ (19) 

The quantum mechanical treatment of these would-be zero-modes restores flavor covariance and 
due to the hedgehog ansatz also rotational covariance. Here we also recognize the advantage of 
the parameterization (^, as it avoids the time-dependent rotation of the flavor variant vacuum 
configuration, ^ = 1. A slightly different parameterization for the flavor rotating meson fields 
is discussed in ref [0. 

We now have to evaluate the action for the field configuration (p!9|). This is most conveniently 
achieved by transforming to the flavor rotating frame q' = A{t)q [Q. This transformation 
induces two additional pieces in the single particle Dirac Hamiltonian [Q 

h = ho + Kot + /isB (20) 

which are due to the time dependence of the collective rotations 

K,, = -tA{t)^ A{t) = J2^'^^- (21) 

a=l ^ 



and the flavor symmetry breaking parameterized by the different constituent quark masses 

/isB = r(3 \ D^^\i + D^^K + (/^88 - i)A8)rt . ( 



22) 



For convenience the angular velocities Vta and the chiral rotation T = cos ® — 'J75T ■ f sin ® 
have been introduced. Furthermore the various isospin invariant components in /isb have been 



Table 2: The fermion contribution to the inertia parameter in the quantization rule (P^l). The 
moments of inertia c? and 0^ are given in GeV~^ while the coefficients a\ and (5\ are dimension- 
less. The up-quark constituent mass m is in MeV. Note that when diagonalizing the collective 
Hamiltonian a piece (~ 1/GeV) originating from induced kaon fields is augmented to the strange 
moment of inertia (? [ffl]. 



m 


2 


"vac 


"L 


P'v 


0^ 


PL 


"iv 


O^lvac 


aitot 


PlV 


/^Ivac 


/^Itot 


350 


7.20 


1.13 


8.33 


1.65 


0.52 


2.17 


-0.925 


-0.004 


-0.929 


-0.203 


-0.003 


-0.206 


400 


4.55 


1.26 


5.81 


1.33 


0.52 


1.85 


-0.401 


-0.005 


-0.405 


-0.121 


-0.003 


-0.124 


450 


3.46 


1.33 


4.79 


1.14 


0.50 


1.64 


-0.232 


-0.006 


-0.238 


-0.081 


-0.004 


-0.085 



made explicitH. The modifications of the single particle Dirac Hamiltonian will be treated 
perturbatively since they are of sub-leading order either in 1/Nc or in {m — 'ms)/m. Details of 
this calculation, which results in a Hamiltonian for the collective coordinates, may be traced 
from the literature |^. The canonical momenta Ra appearing in this Hamiltonian are related 
to the angular velocities 



R„ 



( -{a^na + aiDsa) 

-iP^Qa + PlDsa) 
NcB 
2^3 ' 



a= 
a= 
a= 



=1,2,3 
=4,.. ,7 



(23) 



The inertia parameters a^, ..,/3i are functionals of the soliton profile, again for explicit expres- 
sions we refer to the literature 0. Furthermore Dab = (l/2)tr(AaAA(,yl"l') denote the adjoint 
representation of the collective rotations. They appear in eq ( p3D as a consequence of flavor 
symmetry breaking, i.e. ai and Pi are due to hsB and hence proportional to the mass differ- 
ence m — rris. In table (^ we list the numerical values for these inertia parameters. In the 
process of quantization these canonical momenta are identified with the right generators of 
SU{3). The resulting Hamiltonian H{A, R) whose eigenstates are identified with the low-lying 
I and I baryons, also contains symmetry breaking terms. As a result of exactly diagonalizing 
H{A, R) the ^ baryons cease to be pure octet states, rather they contain sizable admixture 
of higher dimensional SU{3) representations. Similarly the | baryons are no pure decuplet 
states ||2^. A major effect of these admxitures is to reduce the strange quark contributions to 
baryon properties, which are computed using the exact eigenfunctions of H{A,R). This can 
easily be understood because due to the inclusion of symmetry breaking it is less probable for 
the heavier strange quarks to be excited. This effect will also become apparent in our results 
for the structure functions. Finally we should mention that the constraint for Rs restricts the 
permissible SU{3) representations to those containing only states with half-integer spins, i.e. 



fermions |36|. 

2b. The Hadronic Tensor in the NJL Chiral Soliton Model 

•^Throughout this paper we adopt the summation conventions i, j, k, ... = 1, 2, 3 and a, /3, . 



,7. 



The hadronic tensor which enters the cross-section of deep-inelastic scattering (DIS) 

W;l{P,q-S) = ^jd'xe^'^<{P,S\\j;{i),j'Jm\P,S) (24) 

represents the defining quantity for the structure functions. They are certain combinations (see 
below) of the form factors appearing in the Lorentz-covariant decomposition (omitting parity 
violating contributions) 

W;l{P,q;S) = (-g^, + ''j^\Mt,W,{x,Q-') 





+^e^.A.V^ [9i{x,Q')+g2{x,Q^)\ S'^ - ^P''92ix,Q') . (25) 
P-q y J q-P J 

In eqs (p^ ) and ( [25| ) the nucleon state is characterized by its momentum P and spin-orientation 
S. In the nucleon rest frame they are P'^ = {Mn,0) and S*^ = (0, s). The Lorentz in- 
variant quantities are the momentum transfer squared, Q^ = —q^ and the Bjorken variable 
X = Q^/2P ■ q where q denotes the momentum of the exchanged boson. This boson couples 
to the hadronic current J° with a being a suitable fiavor componenta. Here it is important to 
note that this hadronic current is completely determined by the symmetry currents of QCD. 
Hence the hadronic current can unambiguously be computed in any model which has the same 
symmetries as QCD. This model current then completely saturates the hadronic tensor and at 
this point there is no need to identify the degrees of freedom of QCD within the model. 

In the NJL model the interaction among the quarks does not contain any derivatives as can 
be seen from the defining equation (|l]). Hence the currents are formally identical to those of a 
free Dirac theory, i.e. J" = q'yfi{'j5)t°'q with f^ being a fiavor generator. We will concentrate 
on the leading twist pieces of the structure functions. These are obtained from the Bjorken 
limit — g^, 2P ■ q ^ oo with x fixedcl. In this limit the computation of the hadronic tensor 
facilitates considerably since we may adopt the free correlation function for the highly off-shell 
intermediate quark propagating between the two boson insertions: 

g(O7M^"^tVg(0) - g(0)7.t^^t'^7^g(O|P, S) . (26) 

Apparently the hadronic tensor acquires contributions from intermediate quarks propagating 
in forward (^ — > 0) and backward (0 -^ ^) directions. In order to apply this formalism to the 
soliton configuration, which breaks translational invariance, we still need to generate states of 
good linear momentum. This is essentially done in the same way as the fiavor-spin projection 



X PS 



^We omit the flavor indices of the form factors. 

^That is, the higher order contributions in 1/Q^ are dropped. The logarithmic dependence on Q^ originating 
from perturbative gluon emission to the nucleon will be discussed in section 5. 



discussed previously by introducing a collective coordinate R describing the position of the 
classical soliton in coordinate space. As the nucleons in the bra and ket states have the same 
momenta this boils down to an additional integral in coordinates space yielding []T3 

i2ny 



Wf^{P, g; S) = —^ / dt £^i rf^G / d^k sign(A;o) 6{k'^)e'^'^''+^''^'-'^'^+^^<^'~^^^ 



x(iV|g(ei,t)7;.t'^^tVg(6,0) -g(6,0)7.t^|^t'^7^g(ei,t)|iV). (27) 



A^) is defined via (R 



Here, 



P, S) = \/2E exp{iP ■ R) N). Although this now is a simplified 
version of the hadronic tensor it is not well suited for the ongoing calculation in the bosonized 
NJL model. The main reason is that the matrix element to be computed is that of a bilocal 
and bilinear quark operator which appears as an ordinary product rather than a time-ordered 
product. Only the latter would be suitable to undergo the bosonization procedure which in 
turn could provide a rigorous regularization of the structure functions consistent with the me- 
son sector of the model as well as other soliton calculations. One possible way to achieve that 
goal would be to start from the forward Compton amplitude whose absorptive part is identical 
to the hadronic tensor. This amplitude is calculated as the matrix element of the time-ordered 
product of the electromagnetic currents and can hence straightforwardly be implemented in 
the functional bosonization thereby implementing a unique regularization. Although first re- 



sults [^ in this direction are promisingSl we will presently avoid the ambiguities connected 
with regularizing structure functions by employing the valence quark approximation. This 
approximation consists of substituting the valence quark wave-function 

g(e, t) = Mi, t) = e-v*A(t) UyiO + E ^^(0 ^^\^-^ + ^'^\^A (28) 



in the above expression for the hadronic tensor (|2^) . The main arguments to consider this 
valence quark approximation to be reliable are the facts that the valence quark contribution 
does not undergo regularization in any event and the observation that this level dominates 
the results for the static properties. The axial properties are actually almost saturated by the 
valence quark contribution 0. It is important to note that this valence level should not be 
confused with the valence quark in the parton model description of the structure functions. Here 
it rather refers to the distinct quark state which is strongly bound in the soliton background. 
Note that both the rotational and symmetry breaking corrections have been included. This, 
of course, is necessary in order to disentangle various fiavor contributions. Consistent 1/Nc 
counting makes it necessary to consider the time- dependence of A[^. This additional bilocality 



in the time coordinates is treated in an 1/Nc expansion employing the quantization rules (|2 

dA{t) 



A{t) = A{to) + {t-to)- 



dt 




^ Since the incomplete gamma-functions appearing in the proper-time regularization have various unphysical 
poles in the complex plane it is more appropriate to employ the Pauli-Villars regularization when identifying 
the cuts of the Compton amplitude [pq. In contrast to earlier two flavor calculations of the structure functions, 
which only included a single subtraction pT] , E8l p9| , the computation of the Compton amplitude allows one to 
consistently incorporate the two subtractions demanded by the quadratically divergent gap-equation. 

10 



= A(to) + ^(t-to)^(to)E^"^a(to) + ..., (29) 

^ a=l 

where the elhpsis refer to higher order contributions in the angular velocities which according 
to eq (p3[) correspond to higher orders in 1/Nc and/or flavor symmetry breaking. At this point 
we recognize that the computation of nucleon matrix elements associated with the expansion 
([29| ) will also demand the matrix element of fig. This component of the angular velocity is not 
provided by the quantization prescription (ESf). We obtain the relevant information from 



f^s = -iti (A^A(t)U(t)) = tr (\^A{t)[H, A{t)]) , (30) 

with H = H{A, R) being the Hamilton operator for the collective coordinates, its eigenstates 
are the baryons, cf. the discussion after eq (P3|). Employing the defining equation for right 
generators, i.e. [A,Ra] = AXa/2 yields 

(fig) = ^ ^. (31) 

It is straightforward to verify that the application of the same procedure to r2i, .., fiy is consistent 
with the quantization rules (p3D. 

3. Unpolarized Nucleon Structure Functions 

We have now collected all ingredients necessary to apply the valence quark approximation 
to the computation of the hadronic tensor in the NJL chiral soliton model. In a first step we 
project the unpolarized structure function fi{x) from W^,^ by contracting it appropriately with 

^'}>yr^^/^n whore ,= (l + |^)". (32) 

In the Bjorken limit, which in the nucleon rest frame P^ = {Mn, )^ is implemented via 

Qo = \q\ — xMm with |g| — * oo , (33) 

the leading twist pieces of the structure functions are extracted. In this limit this projector 



becomes as simple as —g^^/2 when contracted with the hadronic tensor in the form of eq (pTl) . 



In this limit the evaluation of the contribution due to the bilocality in the time-coordinate 



([29| ) also becomes feasible. The linear t-dependence is treated as a derivative with respect to 



g°. The g° dependence of the structure functions occurs via the Bjorken variable x. Hence 
this derivative is replaced by ■^-^. Accordingly it is convenient to separate the unpolarized 
structure function /i into two contributions: 

/i(x)=/i(a;)+AA(a;) (34) 

with A/i(x) containing the piece due to bilocality in the collective coordinates (123). By def- 



inition this contribution is already linear in the angular velocities, and hence of sub-leading 

11 



order. Consequently the structure function A/i(x) is completely given in terms of the classical 
fields, i.e. by omitting the cranking correction in eq (|28|). 

It is henceforth useful to define the following Fourier-transformations 



and 



lpy{p) 



^vip) 






hsB\y) 



exp (i^-p) 



(35) 



(36) 



_ ^ — # — * 

According to eq (pTf ) we put p = k + g together with the on-shell condition k = \k\ = qoT ^y 
pertinent for the forward and backward contributions to the hadronic tensor. As in the Bjorken 
limit go tends to infinity and qo — k remains finite, we may change the integration measure to 



k'^dVtr = pdpdcj) 



(37) 



with p = \p\ and (p being the azimuth angle between g and p. The lower limit of the p-integration 
is assumed when k and g are anti-parallel, i.e. p™™ = \M]sfX^ey\ =: MAr|a;:[:|. Finally we make 
use of the single quark wave-function containing only finite momentum modes of the order 
p K, e\. Hence the integral receives non- vanishing contributions only when k - 
With these preliminaries the unpolarized structure function /i becomes 



-q 



15, 26 



Jl 



Nc^(N 



TT 



Mn\x-\ 



pdpd(t)^\{p-){l - a ■ q)A^O^''AipY{p_ 



Mn\x+\ 



pdpd<Pij\j{p+){l + a ■ g)AtC«^A^v(p^ 



N 



AA" 



ah 



rab rab 

Jl- ~ Jl,+ 

Nc d 



(3J 



271 dx 



N 



Mm\x-\ 



pdpd(l)¥y{p.){l - a ■ g)AtO^MfigA^§v(p- 



+ 



Mn\x+\ 



pdpd(l)¥y{p+){l + a ■ q)ngX3A^O''^A^Y{p^ 



N 



A/- + AA 



ab 



(39) 



where we have also used that the valence quark wave-function carries positive parity, i.e. 
V'v(~p) = 7o^v(p)- The momentum p± denotes the integrated momentum with the polar 
angle fixed, cos0^ = Mnx±/p. The flavor structure of the integrand is now comprised in the 
matrices O"'^ = f^t^ which are diagonal for all processes in which we are interested 



O'^'' = .^ol + '^sAs + Ms 



Js'^s 



(40) 



These matrices are given by the physical process under consideration as this specifies the flavor 
matrices f^ and t^ in the hadronic tensors. Since we will ignore Cabbibo mixing, the two sets 
(5o, ^3, 5s) and {5'q, 5'-^, 5g) differ at most by the sign of 5-^ and 5'^. 
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We remark the change of sign in the backward moving quark contribution when going 
from /i to A/i. This stems from the different expansion ( pOf ) caused by exchanging the time 
arguments of the collective coordinates. 

According to the quantization rule (|23| ) the collective coordinates and their time derivatives 
Qa are elevated to non-commuting operators. Guided by symmetry requirements such as 
hermiticity, PCAC and G-parity we adopt the sjniimetric ordering whenever ambiguities occurQ, 

i.e. Dab^c -^ {-Dafe,fic}/2. 

For a consistent formulation of the valence quark approximation we consider the sum rules 
for isospin and strangeness of the nucleon. The former is known as the Adler sum rule and 
refers to (anti)-neutrino scattering off the nucleon. For these processes we require f^ = {t^Y = 
(A^ ± ^A^. According to eq ( ^0|) the subtraction of neutrino and anti-neutrino scattering 



processei 



s 



dx ^^p^ _ ^^^r^^ ^^^^ 



X 



corresponds to S'^ = —5'^ = 1 with all other 5J^'s identical to zero. It is then straightforward to 
verify that 

Jdx(f\^-f\+) = {N\2h\N) (42) 

y'rfa;(A/i,_ + A/i,+) = (43) 

with /a = DsiRi + D^aRa + DssRs- However, for eq (^) to be valid, we have to restrict the 
inertial parameters appearing in the quantization rule (p3|) to their dominating valence quark 
contribution. Of course, this is expected. The differences between using the total values of 
these inertial parameters and these valence quark contributions can be employed to estimate 
the uncertainties of the valence quark approximation. Within the same restriction also the 
strangeness sum rule for t"- = t'^ = {1 — -\/3A^)/3 holds as one easily verifies that (for Nc = 3) 

j dx (i\,^ + h,+) = ^(A^|^|iV) (44) 

yda;(A/i,_-A/i,+) = (45) 

with S = {2\^/Nc){Ls — Rg)- An interesting sum rule for the unpolarized structure functions 



is the so-called Gottfried sum rule |^ which addresses electron-nucleon scattering, 

dx 



So = j-^{f?-fT)- (46) 



^While G-parity invariance can be maintained for other orderings when going beyond the valence quark 
approximation this is not the case for PCAC pO| |. 

^Here we do not specify the boundaries of the integrals because they will be subject to the boost into the 
infinite momentum frame, c/. eq (p9|). 
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The structure function /2 is related to the form factor W2 in the defining equation ( pSj) by 
/2 = —q'^W2/2MffX. In the Bjorken hmit the Callan Gross relation /2 = 2x/i is identically 
fulfilled. Using the electromagnetic combination 



O"" = O 



ba 



(A3/2 + A8/2v^)2=:0^ 



we therefore find 



Sg = 2^ f dx(sG{x) + Asci 
Stt J ^ 



X 



(47) 



(4J 



with 



sg{x) 



Mn(N 



Mn\x-\ 



pdpd(j)ip\j{p^)(l - a ■ q)A^X3Ailjy{p. 



Mn\x+\ 



pdpd(l)^l{p+){l + a ■ q)A^\3A^Y{p^ 



N 



and 



Asg(x) 



2dx\ 



Mjv|a;_| 



pdpd(j)¥y{p_){l - a ■ g)AU3AA^fig^v(P- 



+ / pdpdc^¥y{p+){l + a ■ q)X^ngA^X3A^Y{p^ 

'Mn\x+\ 



N 



In the parton model language the assumption of a sea quark distribution which is identical for 
all flavors yields the historic value Sq = 1/3. Recent experiments by the NMC group indicate 
a sizable deviation So = 0.235 ± 0.026 |^. This has often been interpreted as an isospin 
violation. However, this is not the case because a non-vanishing lower component of the quark 
wave-function ipy is sufficient to provide the desired deviation from the historic value even for 
identical up and down quark masses pB| , p5[ . Since the quarks experience a mesonic background 
field this effect may be attributed to pion clouds |^ . 



4. Polarized Structure Functions 

The projector suitable to extract the polarized twist-2 structure function gi from the 



hadronic tensor reads 44 



91 



l(2p.qxS, + ^L-sy 



P-q 



with I 

In the Bjorken limit 
multiplication of 



-4M, 



N 



P-q 



2P-qx-{q- S) 



2 1 



and s II q. 



(49) 
(50) 



its contraction with the hadronic tensor (P7|) is equivalent to the 



AMI/ 



2M. 



N 



q-S 



(51) 
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Taking into account that for gi the nucleon spin s should be oriented along the photon mo- 
mentum g which is accomplished by the choice 

q = s = e^ (52) 

we obtain for this polarized structure function by repeating the computation of the previous 
section 

]\/f /"OO 

TT \ JAfjv|x+| 



vr 



NJ 



= ~9t+-~9t- (53) 

^gf = -^-^{N,s = e, / pdpd<p¥y{p+){l + a3h5X'QgA^O''''A^y{p^) 

2n ax ^ Jmn\x+\ 

+ / pdpd(j)^Up^){l - a,h,A^O'"'AngX3^y{p_) N, s = e,) 

= Agf^ + Agl'_. (54) 

Again we have separated the contributions arising from time dependence of the collective co- 
ordinates (^91) 

gf = ~gf + Agf . (55) 

Of course, the flavor indices have to be chosen according to the considered process. We would 
like to remark that similar to the case of the unpolarized structure functions the decomposition 
([55| ) does not directly correspond to an expansion in 1/Nc and flavor symmetry breaking 
because the valence quark wave-function ip^ contains sub-leading orders as well, c/. eq (pHI). 

For the case of the electro-magnetic interaction O = O^'"^', cf. eq (|47|), the sum rules 
discussed below are straightforwardly verified in the valence quark approximation when the 
previous restrictions to the inertia parameters as well as the static properties are applied. The 
Bjorken sum rule |^ is obtained from the difference between the first moments of the proton 



and neutron polarized structure functions 

T?-T^r^ = Jdx{g['^-g^-^)=lgA (56) 

with gA being the axial (isovector) charge of the nucleon measured in neutron beta-decay. In 
flavor SU(3) the isosinglet charge contains octet as well as singlet pieces. The sum rule for the 
sum of the first moments is obtained as 

r(^) + r(") = ldx {g['^ + ^i")) = ^ ( A« + Ad) + 1 As (57) 

where Ag = {N,s = ez\q'y3'y5q\N,s = e^ refers to the nucleon matrix element of the axial 
current of flavor q. In this notation we have gj^ = Au — Ad. 
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The flavor singlet combination 

Au + Ad + As = J: (58) 

corresponds to twice that part of the nucleon spin which is carried by the quarks, i.e. S = 2Sq. 
Its study has lead to the proton spin puzzle as the EMC [||], SMC and SLAC |8|, H 



measurements combined with the assumption of flavor symmetry and data from semileptonic 
hyperon decays indicated that S was unexpectedly small |ll|]. In addition that analysis yielded 



a large polarization of the strange sea in the nucleon, i.e. a large As/ Ad. However, it was also 
observed that this ratio is quite sensitive to the assumption of SU(3) flavor symmetry, while 
the smallness of S remains almost unchanged once this assumption is waived ED] . Our results 



for Aq have already been presented in [E^. In this context it should be noted that (as long as 



^3 = ^3) the bilocal corrections do not alter Ag, i.e. J dx{Agf\ + Agf'_) = 0. 

5. Projection and DGLAP-Evolution 

The structure functions of the previous sections are calculated in the nucleon {i.e. soliton) 
rest frame (RF). In the introduction we have already argued that we should transform these 
structure functions to the infinite momentum frame (IMF) in order to mitigate the effects 
originating from omitting the dynamical responses of the soliton on the transferred momentum. 
This boost causes a Lorentz contraction [|^, ^ 

/imf(x) = ^Y3^/Rf( - In (1 - ^)) (59) 

which constrains the structure functions to the interval x G [0, 1] hence providing proper sup- 
port. The feature that the structure functions vanish for a; < 1 is not only demanded by Lorentz 
covariance but also mandatory in order to include the logarithmic corrections in Q^ employ- 
ing the DGLAP evolution program |5l|, ^ of perturbative QCD which otherwise contained 



ill-defined singularities. These corrections are incorporated by integrating differential equa- 
tions (c/. eqs (|60|) -(|62|)) for the nucleon structure functions with respect to t = ^(Q^/Aqcj-,). 
The assumption underlying the application of the QCD evolution program to the structure 
functions obtained in the present model (actually any model) is that the model approximates 
QCD at a low scale, Qq. Hence at this scale the structure functions computed in the model 
and those of QCD should approximately be equal. This low scale Qo represents the initial 
boundary to = 1ii(Qo/^qcd) ^^ integrating the differential equations. It should be stressed 
that the lower scale Qg represents a new parameter being intrinsic to the model. We will fix 
Qq by demanding a best fit to the experimental data for the unpolarized structure function at 
the upper boundary of the t-integration. To be specific we will consider the linear combination 
entering the Gottfried sum rule (03). As the DGLAP procedure relies on perturbative QCD 



we should not take too small a value for Qq in order to stay in its range of validity. Of course, 
the upper boundary of the t-integration is set by the experiment to which the model results 
are compared. 



^Hence the structure functions are (reference) frame-dependent quantities. 
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Before actually integrating the DGLAP equations the flavor singlet (s) and non-singlet (ns) 
pieces of the structure functions have to be disentangled as the former may also contain gluon 
contributions. In the notation of eq (HO) the singlet component is denoted by 6q while both ^3 



and 6s are non-singlet. Formally the differential equations to be integrated read 



rf/(-)(x,t) _ «QCD(t) f^dy 



X 



dt 




dt 27T Jx y \y 

df^'\x,t) 



-P,A- r'\y,t), (60) 



dg{x,t) ^ aQCDJt) [^ <^\ 
dt 27r Jx y \ 




(61) 

P,A-\g{y.t) + P,A-\r{y,t)\ . (62) 



Here the quantity / symbolically represents both the unpolarized and polarized structure func 
tions /i and gi, respectively. Note that the splitting functions Pij which are listed in ref pS 



differ in these two cases. At this point we have to make a further assumption as the model 
does not give any information about the gluon content g{x, to) in the nucleon which arises from 
quarks radiating and absorbing gluons and become visible only as the nucleon is probed with 
higher photon momentall2l. We hence assume that at Qg ^^^ gluon distribution function for both 
the unpolarized and polarized structure functions vanish. This is not unmotivated as low-scale 
parameterizations of the nucleon distribution functions indicate that the gluon component is 



indeed small |^6[ ^. Hence the initial values of the integration are the model predictions of 
the nucleon structure functions together with g{x,to) = 0. Note that this assumption is irrel- 
evant when determining to from the unpolarized structure functions which enter the Gottfried 
sum rule (^) because the gluon distribution does not contribute to this isovector combination. 
Having carried out the DGLAP integration to the scale Q^, the singlet and non-singlet pieces 
are recombined to the physical flavor combinations. As a consequence of the evolution (|6^) we 
find non-vanishing gluon contributions to the singlet structure function at any Q'^ > Q^. 

In the numerical treatment we constrain the evolution equations to their leading order 
components in the QCD coupling aQCD(t) = 47r/[(ll — 2A^//3) ln(Q^/AQci-,)] not only for 
simplicityt^ but in particular because a next to leading order calculation requires the identifica- 
tion of valence and sea quark distributions at the scale Q^. This cannot be done in the present 
model as we merely identify the symmetries of QCD rather than the quark degrees of freedom. 
Note again that the valence quark level ( pSf ) is that of the NJL chiral soliton model. It is 
hence a constituent type quark degree of freedom and should not be confused with the valence 
distribution of the parton model. Actually the gap-equation (0) causes the constituent quarks 
to contain not only current quarks but also virtual pairs of current quarks and anti-quarks. 



^^However, gluonic correlations are already contained in the model which can be concluded from the non- 
vanishing twist -3 structure function §2 p6[ . The QCD equations of motion can be used to show that the 
moments of §2 are matrix elements of the gluon field strength [ p4[ . For a review on spin structure functions see 
ref ||. 

^^From bag model calculations it is known that merely the variation of the model scale Qg is the major effect 
of taking the next to leading order correction into account p3] . 
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Table 3: In this table we present our results for the Gottfried sum rule (0). The entries "val" 
and "val+vac" indicate whether only the valence quark pieces or the total values for the inertia 
parameters «^,..,/3i are substituted in the quantization rules (^31). See also table |1|. 





Sg 


ASg 


Sg + A5g 


m (MeV) 


val val+vac 


val val+vac 


val val+vac 


400 
450 


0.29 0.22 
0.27 0.20 


-0.07 -0.05 
-0.11 -0.08 


0.22 0.17 
0.16 0.12 



6. Discussion of Numerical Results 

We are now prepared to compare our model predictions to the experimental data obtained 
for DIS at CERN and SLAG. 

Within this model a reasonable description of the static properties of baryons is achieved 
for constituent quark masses in the range m = 400 . . . 450MeV [Q . We will henceforth consider 
the two limiting values for the only free parameter in the baryon sector. 

We set the stage by discussing our results for the Gottfried sum rule (^6|) in table ||. The 
value of this integral remains unchanged in the process of transforming to the IMF as well as 
the subsequent leading order DGLAP evolution. In presenting the results we have disentangled 
the two pieces associated to the integrand in eq (H). Apparently the bilocal corrections provide 
a sizable negative contribution; a result which was already observed in the two flavor model 
39| . In table |^ the numerical results for two different calculations are shown. Those labeled by 



'val' originate from substituting only the valence quark piece of the inertia parameters when 

eliminating the angular velocities in favor of SU(3) operators (|23|) and matrix elements (0). In 

this case consistency of the valence quark approximation with the Adler sum rule is maintained 

26|, see also section 3. The entry 'val+vac' refers to the inclusion of the usually significantly 



smaller vacuum pieces in the inertia parameters. For nucleon properties which are obtained 
from matrix elements of vector currents the valence quark commonly contributes to about 80% 
to the static properties which in the case of vector charges are evidently related to the inertia 
parameters. For axial properties the valence level more or less saturates the static properties 
1^, |59|. For those reasons we consider the valence quark approximation to be most reliable when 
restricting to the 'val' prescription in the case of the unpolarized structure functions but using 
the total inertia parameters when later discussing the polarized structure functions. After these 
general remarks we return to the Gottfried sum rule and observe that for m = 400MeV a fair 
agreement is obtained with the empirical value 0.235 ± 0.026 ^2|. Our result for m = 450MeV 
is apparently on the low side. This may be due to the valence quark approximation getting 
worse for larger coupling because a stronger binding of the valence quark is joined by a more 
pronounced polarization of the vacuum. 

We proceed by comparing the integrand of the Gottfried sum rule, sg + As^ in eq (|48|) , with 
the empirical data, which we assume to obey the Callan Gross relation, /2 = 2x/i. For this 
comparison we apply the DGLAP evolution from the model scale Qq to the scale Q^ = AGeV^ 
at which the data are available. At Qq the integrand of the DGLAP equations is taken to be 



Figure 1: Numerical results for the unpolarized structure functions entering the Gottfried sum rule 
for m = 400MeV (left panel) and m = 450MeV (right panel). Here the label RF indicates the 
structure functions as calculated from the expressions in the appendix, the transformation to the 
infinite momentum frame ( |59D is indicated by IMF and finally the structure functions resulting from 
the leading order (LO) DGLAP evolution are represented by the solid line. The empirical data are 
those for (/I'' — /™) /2a; |0. Only the valence part of the moments of inertia enters the quantization 
rule. 
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the combination sq + As^ transformed to the IMFEJ. By assimilating to the overall shape of 
the data we determine the model scale to be Qq ~ 0.5GeV , as shown in figure |l|. A smaller 
value for Qq might even improve on the agreement as it lowers the structure functions at large 
and moderate x by simultaneously increasing it at small x. However, for too small a Qq one 
would have to go beyond leading order in the splitting functions. As already discussed above 
this would require additional assumptions when identifying the model structure functions with 
those of QCD. From figure |I] we also note that at x ~ 0.7 the combination entering the Gottfried 
sum rule becomes slightly negative. This negative contribution is completely due to the bilocal 
corrections Asq- 

Now we turn to an interesting quantity only accessible in the three flavor extension of the 
model: the contribution of strange degrees of freedom to the unpolarized structure function, 
fl{x). Of course, a major motivation for the present study is the investigation of quantities 
hke fl{x). This structure function is obtained by setting C"^ = O'"^ = (1 — y/SXs)/^- Note that 
the integral over /f (x) is not given by the strangeness charge discussed in eq ( ^41) as for these 
two quantities the pieces stemming from the backward moving quarks have opposite signs. Our 
results for /f (x) are shown in figure ^j. 

We observe that the structure function /f exhibits a pronounced maximum at x ~ 0.5 after 
boosting the rest frame structure function to the infinite momentum frame. Naively one expects 
the momentum distribution at the low energy scale to be peaked around x = 1/3. However, 



^^Although the Callan Gross relation holds in our model for the structure functions at the model scale Qq 
it does not necessarily hold after projection and evolution. Hence we consider it appropriate to compare our 
results with the data obtained by applying the Callan Gross relation to the available data for /2. 
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Figure 2: Model predictions for the unpolarized structure functions of the strange projector for two 
different values of the constituent up quark mass. Only the valence part of the moments of inertia 
enters the quantization rule. 
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the constituent mass of strange quark is about 50% larger than that of the up or down quark, 
cf. table 0, which explains the shift of the maximum to a larger value in x. We also note 
that this maximum increases with the constituent mass m. This feature can be understood by 
recalling that this parameter actually measures the coupling of the non-strange quarks to the 
soliton with the soliton then exciting the strange quarks. The DGLAP evolution only leads to 
a moderate smearing of the maximum. 

At lower x the strangeness projected unpolarized structure function becomes negative which 
makes a probability interpretation in the spirit of the parton model impossible. This short- 
coming is most likely linked to the valence quark approximation which not only omits contribu- 
tions from the polarized vacuum (which still miss the formulation of a consistent regularization) 



as well as certain disconnected diagrams |T5[ which in the background of the soliton are not 
necessarily disconnected. Apparently the property of a negative piece in fl{x) gets mitigated 
at higher scales as the DGLAP evolution increases this structure function at small x, it even 
might diverge as x ^ 0. We would like to remark that this divergence is already present when 



the bilocal corrections A/i are omitted pO 



As already pointed out in the introduction the polarized structure functions are of particular 
interest. Here we present the results for both electron-proton as well as electron-neutron 
scattering in figures ^ and |^, respectively. We see that after projection and DGLAP evolution 
the results for electron-proton scattering reproduce the gross features of the empirical data 
49| for both values of the constituent quark mass m = 400MeV and 450MeV. However, at 
moderate x the predictions slightly overestimate the data while at larger x they drop off too 
quickly. This might be linked to the boost into the IMF which strongly squeezes the structure 
functions as the comparison between the RF and IMF structure functions in figure ^ indicates. 
This is especially the case for smaller constituent quark masses as for those the RF structure 



functions are already moderately localized. Other projection procedures [31] may have a weaker 
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Figure 3: Numerical results for the polarized structure functions for electron-proton scattering for 
two different values of the constituent up-quark mass. The experimental values are taken from the 



SLAC E143 experiment ||4^. The sum of the valence and the vacuum parts of the inertia parameters 
enter the quantization rule. 
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effect. We would also like to mention that the bilocal corrections Agi{x) tend to enhance this 
pattern. 

By taking appropriate matrix elements in eqs ( [^ ) and ( p^ we straightforwardly obtain the 
polarized structure functions for electron-neutron scattering. When comparing these results to 



Figure 4: Same as figure |]for electron-neutron scattering. 
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the data of the E143 experiment [^ in figure ^ we observe that the model properly reproduces 
the small trough in the experimental data at x ~ 0.5. On the other hand our model calculation 
does not share the tendency of most of the data to become sizable and negative at smaller x. 
To us it seems that the errors on the data are still too large to draw very stringent conclusions 
in this respect. 

Finally we turn to the strange quark contribution of the polarized structure function. Again 
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this is defined by setting 0°-'' = 0^°- = (1 — v^A8)/3. Our results are shown in figure ||. These 
have aheady been presented previously [^, however, with the bilocal correction /\gl{x), which 
is sub-leading in 1/Nc, omitted. The general shape of gl{x) + Agl{x) is not altered by this 



Figure 5: Same as figure ^ for the polarized structure functions. Here both the valence and the 
vacuum part of the inertia parameters enter the expressions for the angular velocity ( ^ , ^11) . 
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additional piece. In particular the significant variation with the constituent quark mass m is 
maintained. We observe from figure |^ that this structure function has the strong tendency 
to become more negative at larger x with increasing m. In any event, we should emphasize 
that the magnitude of the strangeness contribution to gi{x) is quite small for all permissible 
parameters and the conclusion of a more or less vanishing polarization of the strange quarks in 
the nucleon is conceivable. The conclusion that strangeness degrees of freedom yield at most 



Figure 6: Comparison with experiment [^ of the two and three flavor model predictions for the 
polarized structure function gi{x) of electron nucleon scattering. 
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moderate effects can also be drawn from figure M where we compare our present results for the 
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polarized structure functions to those obtained in the two flavor reduction of the same model. 
The isoscalar combination g^^ + gf^ remains almost unchanged by the inclusion of strange 
quarks while there is a slight increase for the isovector combination gl^ — gf. The fact that the 
latter changes at all is not only due to the additional terms (c/. appendix C) when extending 
the model to three flavors but also to the change between SU(2) and SU(3) Clebsch-Gordan 
coefficients arising in the calculation of nucleon matrix elements of collective operatorsli^l. We 



already communicated this result in ref p9[ for the case of the polarized structure function 
gi{x). As can be observed from figure 1 in that reference, the use of fiavor symmetric collective 
wave-functions enhances the difference between the extrema of gi{x) by almost a factor two. 

The smallness of the strange quark contributions to the nucleon structure functions is (at 
least partially) due to fact that we employ a nucleon wave-function obtained by exactly diag- 
onalizing the collective Hamiltonian. This mitigates the role of strange quarks in the nucleon 



T0| . We have checked that using an SU(3) symmetric wave-function yields larger strange quark 



pieces in the nucleon structure functions. 

Figure 7: Gluon contributions to the unpolarized (left panel) and polarized (right panel) structure 
functions as they arise from the DGLAP evolution ( |62D . Note the different scales for the unpolarized 
and polarized contributions. 



Gluon Contribution to f j 

Q,,' = 0.5 GeV", q' = 4.0 GeV' 



m^ = 400MeV 

m° = 450MeV 




Gluon Contribution to g^ 

Q,,' = 0.5 GeV', q' = 3.0 GeV' 



m„ = 400 MeV 
m,, = 450 MeV 



For completeness we also show the gluon contributions to the structure functions at the large 
momentum scale Q^ as they emerge from the DGLAP evolution (|62D in figure ^ We recall that 
as an assumption these contributions were put to zero at the model scale Q^. Apparently 
the gluonic piece in the unpolarized structure function is almost two orders of magnitudes 
larger than in the polarized case. This huge gluonic contribution to the unpolarized structure 
function due to the DGLAP evolution is a well known effect. On a double logarithmic plot 
the unpolarized gluon distribution (times Bjorken x) grows approximately linearly as x ^ 
for large enough Q"^ |^. The dependence of the gluonic pieces on the constituent quark mass 
apparently is negligible. 



^•^When constructing the exact eigenstates of the collective Hamiltonian we employ the full inertia parameters 
when applying the valence quark approximation to the structure functions. 
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7. Conclusion 

We have studied the nucleoli structure functions within the three flavor Nambu-Jona- 
Lasinio chiral soliton model which in particular allowed us to compute the contribution of 
strange quarks to the structure functions. 

Like in QCD, in this model the flavor symmetry is explicitly broken by different current 
quark masses. By the mass generation mechanism of spontaneous breaking of chiral symmetry 
this also induces flavor symmetry breaking for the constituent quark masses. In the baryon 
sector of the model the flavor symmetry breaking effects are explicitly incorporated by using 
an extension of the Yabu-Ando approach to quantize the chiral soliton. This yields a collective 
wave-function for the nucleon with good flavor and spin quantum numbers. In addition a 
projection onto good momentum states is performed to restore translational invariance which 
is broken by the soliton. Having obtained states of good momentum, the structure functions 
computed from a localized object are boosted to the inflnite momentum frame. 

Although for the construction of the chiral soliton the inclusion of the Dirac sea is mandatory, 
in the evaluation of the pertinent matrix elements of the hadronic tensor we use the valence 
quark approximation for the nucleon wave-function. This not only avoids the unsolved problem 
of consistently regularizing the vacuum contribution but should also be sufficient since in the 
chiral NJL soliton model the static properties of the nucleon are dominated or even saturated 
by their valence quark contribution. In this context we should emphasize, however, that in 
our model the valence quarks are bound constituent quarks and thus differ essentially from 
a valence (or sea) quark in the parton model. To minimize the ambiguities when identifying 
model and QCD degrees of freedom at the low energy scale where our model is supposed to 
approximate QCD we only use the leading order DGLAP formalism to evolve the predicted 
structure functions to the energy scale encountered in the corresponding experiments. 

In the numeric evaluation we adjust the parameters of the model to well-known mesonic 
properties leaving only the constituent up quark mass as a free parameter. Note that this also 
determines the constituent strange quark mass. The Gottfried sum rule is fairly well reproduced 
for a constituent quark mass oi m = 400MeV but somewhat underestimated for m = 450MeV. 
This underestimation partially reflects the fact that for large constituent quark masses the 
contribution of the distorted vacuum becomes non-negligible (In this context one should note 
that the constituent quark mass represents the coupling constant between the quarks and the 
chiral fleld). The predicted unpolarized structure functions which enter the Gottfried sum rule 
overestimated the experimental data at large and moderate x but underestimate them at small 
X. A better agreement could probably be obtained by lowering the reference scale of the model 
in the DGLAP evolution. However, then one would have to go beyond the one leading order 
splitting functions. 

A main focus of the present paper has been the contribution of the strange quark to the 
unpolarized and polarized structure functions. The strange quark part to the unpolarized struc- 
ture functions fl{x) exhibits a pronounced peak at x = 0.5. The shift from the naively expected 
peak at a; = 1/3 can be attributed to the larger mass of the strange quark. The DGLAP evo- 
lution leads only to a moderate smearing of the peak. At lower x the strangeness contribution 
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to the unpolarized structure functions is negative which is probably due to the valence quark 
approximation which ignores certain (in the vacuum) disconnected diagrams which in the pres- 
ence of the chiral soliton field become connected. In general, however, the regularized vacuum 
contributions need not restore positivity. The DGLAP evolution substantially increases the 
strangeness contribution at smaller x. 

After projection and DGLAP evolution the calculated polarized structure function for deep 
inelastic electron-proton scattering reproduces the gross features of the experimental data, 
although our results slightly overestimate the data at moderate x while they drop off somewhat 
too quickly at larger x. For electron-neutron scattering our calculated polarized structure 
function properly reproduces the small trough in the experimental data at a; ~ 0.5 but fail 
to become sizable enough at lower x. Unfortunately the error bars on the data seem to be 
too large to allow for very stringent conclusions. Again, choosing a smaller reference scale of 
the model might improve on these deviations from the data. Our numerical results show that 
the strange quark contribution of the polarized structure functions is quite small and in fact 
conceivable with a vanishing polarization of the strange quarks in the nucleon. The comparison 
of the results from the two and three flavor chiral soliton model leads us to the conclusion 
that strange degrees of freedom yield at most moderate effects to the structure functions of the 
nucleon. 
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A. Basis wave-functions 

Although the basis states to diagonalize the quark Dirac Hamiltonian (|I8|) in the back- 



ground of the static soliton may be traced from the literature (c/. ref |^2|) we repeat them 
here for completeness and use in the proceeding appendices. The eigenfunctions of this static 
Hamiltonian are eigenfunctions of the grand spin operator 

where r/2 denotes the isospin operator in the fundamental representation. It is easily verified 
that G commutes with single particle Hamiltonian (|l^). Additionally, the wave-functions are 
eigenfunctions of the parity operator. To describe the coupling of orbital angular momentum 
/, spin j = I ± 1/2 and isospin r = 1/2 to the total grand spin G with projection quantum 
number M we adopt the notation 

yF/\f) = mCM) with \ljGM) = Y: C^JljJ^^ll^ ■ (64) 

The coordinate space wave-functions for the (non-strange) up and down quarks are given by 

' ^^ ~ \ ff'^-^'\r) y^l{a^,,,{f) ) + 1, -/f '+^^)(r) y^'ia-y^i^' ' ^^^^ 
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r 






(66) 



The superscript ± refers to the intrinsic parity Hint which is defined via total parity and 
the grand spin 11 = nint(— l)*^- To discretize the quark-orbits we demand that the upper 
components of these spinors to vanish at the boundary of a spherical cavity with radius DM. It 
should be noted that the classical part of the valence quark wave-function dwells in the G = 
channel with only j = 1/2 being permitted: 

*vW-Uv(r) yZ/.{r))- ^''^ 

As a shorthand notation to deal with the symmetry breaking corrections we introduce here the 
so-called chirally rotated radial wave-functions 

5,®= = ^ cos Y + / sin Y and f^" = / cos -^ - 5- sin ^ (68) 



Here Gc denotes the chiral angle which minimizes the classical energy ([I7|) . Since the classical 
soliton is embedded in the isospin subspace, the strange quarks wave-functions are those of 
a free Dirac field in a spherical cavity. We choose spinors with good angular momentum (c/. 







'i'W'^iJii.knir) {r\ll + \m) \ 

w~iji+i{knir) {r\l + ll + \m) ) 

iw^iji{knir) {r\ll-\m) 
-w-iii-i{knir) {r\l - 1/ - \m) 



with 



Eni = ±Jkli + ml Mnl 



D^I^\ji+,{k^iD) 



'^nl = V 1 + ^s/Enh w^i = sgn{Eni)\/l - ms/E^i . 

Although the diagonalization of the static quark Dirac Hamiltonian operator is most con- 
veniently performed in position space we will need the wave functions in momentum space for 
the computations of the structure functions. The grand spin structure will be preserved under 
the Fourier transformation. Generally we denote by R{p) the Fourier transform of the radial 
function function R{r) keeping in mind that it implements a spherical Bessel function of the 
order of the associated orbital angular momentum, e.g. g'f^'~^'^\p) = Jq r'^drgj^'~^'^\r)JG{pr). 
For the strange quarks the situation is a bit more involved and we find it suitable to define 
Sni{k) as the Fourier transform involving ji{knir) while t\i{k) is obtained from integrating over 
ji+i{knir) and ili{k) contains ji-iiKur). 



"'^'^ Obviously that radius is chosen to be a multiple of the typical soliton extension. 
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B. Valence quark wave-function: 
cranking and symmetry breaking 

Here we make explicit the corrections to the classical valence quark wave-function which 
are due to both the collective rotation /irot and the flavor symmetry breaking part hsB of the 
perturbation in eq (pOf ) (see also eq (pSD). 









V2 



(Dsi - iD 



82 J 



71" 



(Dgl + ^^82) - ^MoD^ 



83 



+ (5Gor^(D88 - 1) 



^ZO^.i 



1 



ev-e„i ^2\^2 

771 — 771 



^^,^(^4 + ^^5) - ^m,i(^6 + ^^7) W^igyiknl) + W„Jy{knl 



'J 2 



V^ 



Sm,^{Dsi + iDss) - S^iiDse + iD^^) w^Cg^^{Ki) + <J^^(A;„0 



"J — 



The summation index yU runs over /i, G and M = —G, .., +G when coupling to a non-strange 
spinor and over n,j = l± 1/2 and rrij = —j, ..+j in the other case. We have used the following 
abbreviations for the overlap integrals: 



Q, = -^^jdTT^gygf'-^ + hf, 



P„ 



ey- 


(^fiJ 


m 


-TUs 


V3(ev - £/.) 


m - 


-nis f 


3(ev 


-e,)J 



(2,-) 



^?^^?-^^^ - /y^^4^'-^^ 



while the momenta kni are the roots oi ji{kniD). 

C. Expressions for the structure functions 

In this appendix we present the expressions entering the structure function calculation using 
the wave-functions of the preceding appendices. 

C. 1 Explicit expressions for /i 

In this part we give the explicit expression for the unpolarized structure function /i in 
terms of the radial wave-functions that have been described in Appendix A. Again, g®" and 
f^" denote the chirally rotated wave functions. 

/i = 5',fl- + ^3/1,- + ^s/iV - 5o/i% - hfl+ - Ssfl+ with 



fU 



TT JMj^\x±\ 

A/f ^ /'f^ ~ ~ 



TT 



a=l 



for / = 3, 



Mjv|x±| 
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N< 



M, 



N 



c- 



TT JMn\x±\ 



pdp ( 9y =F 2/vfi'v cos + /< 



/[ 






+ fv f L''^^T cos e'^igyfJ'+^ + f. 



± 



N, 



c- 



gyg], ' ' + Jvr^ • " ^ cos w (^^vj;^ ' ' + Jvg^, 

pdp I (qw =F 2/v^v cos 0^ + /^ 



n(l'+) 



(iV|(D88-l)|iV) 



Vr JMjv|x±| 



(iViAsliV) 



2^3 



^v^i^'+) 



/v/^+)tcos0±(^v/^+^ + / 



rvfi- 



n(l.+) 



(iV|A8p88-l)|iV) 



+bGi \Q, (b, t 2iV^cos0±) {N\Dunm + 2P^ (P^ t 2N^cosQ^) {N\DuDs^\N) 



1, 



A/'2 
■'^nO 



WnOfi'v(/i;no)SnO + ^i^„0^„o/v(^no)SnO ) (fi'V ± /v COS 



4 ''^ ev-e^ 

+ {'^nO'^no9v{kno)ilo + ^'q fy {kno)ilo) (/v ± ^V COS 0^) | (iV| Aa^a |iV) 
+ 1 («^rlo^^v'(^no)SnO + «^+oW^o/v'(^no)Sno) (^V ± /v COS 0^) 
+ {wnoWno9Y'{kno)il,o + W~o fY^iknoWno) (/v ± OY COS0^' 



x—{m - rris) {N\Di^D8a\N) 



for / = 3, 8 with the summation conventions i = 1,2,3 and a = 4, . . . , 7. In case not exphcitly 
stated, the integration variable p is the argument of gv and /y. We have furthermore used the 
following shorthand notation: 






fvf!'-^ 






;,(2,-) 



in addition to the quantities Q^, P^ and F^ which were already introduced in appendix B. 

C.2 Explicit expressions for A/i 

As usual we split up the structure functions according to the flavor content and the two 
contributions from the forward and backward intermediate quark. In the case of A/i this yields 



A/i = <5oAA% + SsAfl^ + 6sAfl^ + 5^A/o_ + 6',Afl_ + 6'sAfl 



with 



A/?, 



Nc d f°° 



± 



27r dx Jmn\x±\ 



pdp (h{pf + gwipf =F 2fy{p)gy{p) COS0^ 



1 



la^ 12/^2 
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while for / = 3, 8 we obtain 

Nc d r°° 



A/U 



27r dx Jmn\x±\ 
1 



pdp (fvipf + 9Y{pf =F '2fv{p)gv{p) cos 6= 



^^-{N\DuiR^ + a,Ds,)\N) - -^{N\DURa + (3iDsa)\N) 



+- 



(iViAsliV) 



C.3 Explicit expressions for gi 

In this appendix we give the exphcit expression for the polarized structure function gi in 
terms of the radial wave functions that were described in Appendix A. As previously g®" and 
f^" denote the chirally rotated wave functions. 

gi = Sogl+ + 63gl+ + S8gl+ - ^ofi-?, - - ^sfi'i - - ^sfi'i - with 

A/f /'C'O _ _ 



aU 



N> 



c- 



TT 



Mn\x±\ 
8 



sL 



9U 



Nc^{N\J2Dia i^ , ^pdpd(j)^l,{p±){l±a3h'>XaMP±)\N) 



TT 



a=l 



■,t 



for / = 3, 



Mm\x±\ 






Qr 



+Pa 



^ (v^Ar(^) + Ar(^)) t ^ (v^M(^)^ + M(^)^ 



£/i,± 



iV, 



M 



Af 



C 



vr JMjv|z±| 



prfp 



^v/v COS 0=^ T 2^v T (cos^ ©^ - 2 )/v 



{N\n,\N) 

{N\Dss\N)\ 
{N\Dis\N) 



-2 E ^Gor, 



COS 6= 



1 



^v/i^'+) + /v^i^'+)) T -gvg^^^^ T (cos^ e± - ^)/v/(^'+) 



1. . .. 

2' 



x(Ar|A3p88-l)|iV) 



7!^ 



+— Qm 



1 X' 






(A^|A8fi3|A^) 
(Ar|A8l)83|A^) 



4ev-e 



nO r 



^vcose±T-/v(2cos2e±-l)) 

+ 2 (/v COS 0=^ =F ^v) {wj,^gM{kno)SnO + W^W^^o/v(/i;no)SnO 

x(Ar|2rf3a/3Aafi/3|Ar) 
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m — rris 



M? 



nO 



2v^ ev-e^ 



^vcose±T^/v(2cos2e±-i; 



^nO 



X ['WnoWno9v^ikno)ilo - Wj f^^{kno)t 
+ 2 (/VCOS0=^ =F ^v) {wno9Y''{kno)SnO " WnoWnofv" {^no)SnO 

x{N\2d3al3Di^Dst3\N)]\ 



for / = 3, 8. The summation involves i = 1, 2, 3 as well as a, (3 = A, ... ,7. We have used the 
following shorthand notation in order to simplify the presentation 



nU) 



?(i,-) 



avflf' ^ + fv~9^ 



nO'-) 



J = 1,2 



<^^ = M^-^S cos^ 0± - 1) + /v/(^'-)) (l + cos^ 0^) 






2(^v^i^'-) + /v/i^-H2cos2e±-l 



while Q^, P^ and F^ are again those listed already in appendix B. 
C.4 Explicit expressions for /S.gi 

Finally we present the formulae for the bifocal correction Af^i. 

Herein the expression for the singlet is quite simple 



A^U 


i\C ^ 1 


2tt dx Jmn\x±\ 


while for / = 


= 3, 8 we find 


^9U = 


Nc d f°° 

Zir ox Jmn\x±\ 



pdp 



1 



1. ^ 



cos e^Mp)Mp) T ^Mp)' t (cos^ 0=*= - -)Mpf 

^ -^{N\Ds3\N) 



2a2 a2 



cos0=t^v(p)/v(p) T IMp? T (cos^ e± - ^)Mpf 



X 



i(7V|.«fl3„fl.|A'>-(^ + ^-^ 



(A'ICasIA') 



-;)3(-|+5)<-|^=»^-|-)+ii^ 



{N\Ds3Dss\N) 



\^a' 



:{N\DssR3\N) 
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A^U = iTir P^P cose^gy{p)Mp) T 7;Mp? T {cos' e^--)Mpy 

Zn ox JMn\x±\ I I I . 

X - i(iv|d3„.fl,„fl.iiv> - (^ + 6^ - ^) mo^w 

This completes the hst of matrix elements entering our structure function calculation. 
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